
function y = difNSD_MultiObj(w,t,P,X,Z,R,C)

y = zeros(t,1);
PZZ = calculatePHat_fixedW(w,t,X,Z,R,C);
for i = 1: t
  
    PP(:,:) = P(i,:,:);   
    PZZZ(:,:) = PZZ(i,:,:);
        
    for ii=1:R
        for jj=1:C
            y(i) = y(i) + myfunc(ii,jj,PP,PZZZ);
        end
    end
     
end

end



function v = myfunc(ii,jj,P,PZ)

p = P(ii,jj);
pHat = PZ(ii,jj);

if pHat ~=0
    v = (p-pHat)^2/pHat; 
else
    v = 0;
end

end

